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Abstract 

An ideal of a local polynomial ring can be described by calculating a standard 
basis with respect to a local monomial ordering. However the usual standard basis 
algorithms are not numerically stable. A numerically stable approach to describing 
the ideal is by finding the space of dual functionals that annihilate it, which reduces 
the problem to one of linear algebra. There are several known algorithms for finding 
the truncated dual up to any specified degree, which is useful for describing zero- 
dimensional ideals. We present a stopping criterion for positive-dimensional cases 
based on homogenization that guarantees all generators of the initial monomial ideal 
are found. This has applications for calculating Hilbert functions. 

A Grobner basis for a polynomial ideal provides a wealth of computational information, 
for example the dimension of the ideal, its Hilbert function, a way to answer the ideal 
membership question, and more. Computing a Grobner basis is a well understood problem 
at least in the setting of exact computation, for example using the Buchberger algorithm. 
Roughly the same can be said about ideals in a local polynomial ring. In the exact 
setting we can compute a standard basis (the local equivalent of a Grobner basis) using 
variations of Buchberger's algorithm such as those using Mora's Normal Form algorithm. 
A treatment of ideals in local rings and standard basis algorithms can be found in [3] and 

However in many practical situations, using only exact computations becomes infea- 
sible and we are forced to rely on approximate numerical data. For example many large 
systems of polynomials can only be solved in practice with numerical algorithms such as 
homotopy continuation. We may want to investigate the properties the ideal in the local 
ring at some solution point, but this point is only known to us approximately. Although 
we can approximate the point with arbitrarily high precision, the error can never be en- 
tirely eliminated. In this context the usual algorithms for computing a standard basis are 
unsuitable because they are not numerically stable. Even arbitrarily small errors in the 
initial data can produce results that are combinatorially incorrect, for instance incorrect 
values of the Hilbert function. Some of the alternative approaches for computing Hilbert 
functions, such as Janet basis algorithms [I], must also be ruled out because they lack 
numerical stability. Other approaches useful in the exact setting can be found in [TT] . 

In a numerical setting, to compute the information provided by a Grobner basis we 
need to tread carefully because many tools are no longer available to us. One avenue 
developed by Bates, et al. |2] is to find witness points of the various components of 
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the variety in order to compute dimension of the ideal and other information. Another 
approach that can be used in the case of zero-dimensional ideals is computing a border 
basis of the ideal, which has better numerical stability than Grobner basis computations 
[7J. In this paper we will focus on an approach which uses purely local information: 
computing the local dual space of the ideal. The dual space is the vector space of all 
functionals that annihilate every element of the ideal. This idea was first developed in 
the seminal work of Macaulay [10J. The dual space of an ideal can provide much of the 
same information as a standard basis, such as the Hilbert function of the ideal, and a test 
for ideal membership. 

There are several algorithms for computing the dual space of an ideal in a local ring, 
truncated at some degree. One that will be discussed in this paper is the Dayton-Zeng 
algorithm presented in [4j, and another is the Mourrain algorithm presented in [12], al- 
though both are based on the ideas of Macaulay. The numerical advantage to dual space 
algorithms is that they reduce the problem to finding the kernel of a matrix. This can be 
done in a numerically stable way using singular value decomposition (SVD). 

These truncated dual space algorithms provide a way to fully characterize the local 
properties of an ideal when the ideal is zero-dimensional, i.e. the point of interest is an 
isolated solution. In this case the dual space has finite dimension, so truncating at a high 
enough dimension we will find a basis for the whole space. However, when the ideal is 
not zero- dimensional or when the dimension is not known a priori, this strategy will fall 
short. 

Our contribution is a method of finding the truncated dual space up to sufficient de- 
gree to ensure that the important features of the local ideal are found. In particular, 
this means finding an explicit formula for the Hilbert function of the ideal at all values, 
not just the values up to some finite degree. The method presented can also be used 
to recover a standard basis for the ideal, which as far as we know is not possible us- 
ing existing truncated dual space algorithms alone in a numerical setting. Additionally 
these tools can be used to answer the ideal membership test for polynomials up to some 
bounded degree. In this way we can describe the local properties of an ideal numerically, 
using purely local information, for an ideal of any dimension. We have implemented 
this method in the Macaulay2 computer algebra system [5j. Our code can be found at 
http : / /people . math . gatech . edu/~rkrone3/NHcode . html , 

A potential application for this result is for developing numerical algorithms for com- 
puting the primary decomposition of an ideal. Current algorithms for primary decom- 
position use elimination theory which relies on Grobner bases. On the numerical side, 
there are algorithms for decomposing a variety into irreducible components. As discussed 
in [8] , this is done by intersecting with random afline spaces of the correct dimension to 
collect witness points on various components. Then homotopy methods are used to de- 
cide which witness points belong to the same components. This is part way to a primary 
decomposition, since the irreducible components correspond to the minimal associated 
primes. However the remaining obstacle is numerically detecting embedded components 
of the ideal. Given a point of interest on the variety that has been found numerically, 
knowledge of the Hilbert function may help decide whether or not the point sits in an 
embedded component. 

In Section [T] we describe preliminary information, defining a standard basis and the 
Hilbert function of an ideal in a local ring, as well as algorithms to calculate them. In 
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Section [2] we define the dual space of an ideal, and show how the dual space can be used 
to recover information about the ideal. We also describe the Dayton-Zeng and Mourrain 
algorithms here. In Section [3] we show how the homogenization of an ideal I motivates 
an algorithm for finding the truncated dual up to sufficient degree and Sections @] and \5\ 
contain our main result: an algorithm for finding the Hilbert function and standard basis 
for an ideal in a numerically stable way using the dual space. 

1 Preliminaries 

Let fi, . . . , f s be a system of polynomials in ring R = k[x\, . . . , x n ] where k is a field, and 
let J = (fi, . . . , f s ). In practice we can assume k = C because we are interested primarily 
in numerical applications. Suppose the point b G A n (k) is known to be in the zero set of 
these polynomials, but b has been calculated numerically so it may not lie exactly on the 
variety. We would like to characterize the zero set V(J) in a neighborhood of b. 

The proper context for answering these local questions is in the local ring at b. Let Rb 
be the localization of R with respect to the maximal ideal m = (x\ — b\, . . . ,x n — b n ), so 

R b = {f/g\f,geR, g{b)^0}. 

Let I be the extension of J in this local ring / = JR^. Without loss of generality we 
will take b = (0, . . . , 0). This makes calculations simpler, and for b ^ we can translate 
elements of Rb to the origin by substituting each %i with Xi + 6j. Every / G R Q can be 
expressed as a power series which converges in some neighborhood of the origin, so 

/ = J2 c - x ° 

a£N" 

with each c a G k. Here a = (a%, . . . , a n ) is a multi-index and x a denotes the monomial 
x" 1 • ■ ■ x" n . Let \a\ — a% + ■ ■ ■ + a n , i.e. |a| is the degree of x a . 
The ring R Q is equipped with a local order >. 

Definition 1. A local order is a total order on the monomials of a local ring that is 
compatible with multiplication, and has 1 > Xi for all Xi (in contrast to a monomial order 
where 1 < Xi). 

Taking the reverse of any monomial order produces a local order and vice versa. 
We will take the local order to be anti-graded, meaning that it respects the degree of 
the monomials, similar to a graded monomial order. Let in> / denote the lead term of 
/ G Rq. Note that even if / is not polynomial, it still has a well defined lead term when 
considered as a power series. Let in> / denote the initial ideal of /, in> / = (in> f\f G J). 

In an exact setting questions about I could be answered by finding a standard basis, 
which is the local equivalent of a Grobner basis. 

Definition 2. Given a local order > on R , a standard basis G of ideal / is a finite set 
G={gi,...,g r }Cl with (in> g 1: . . . , in> g r ) = in> I. 
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The simplest algorithm for calculating a standard basis is the Buchberger algorithm 
(as with a Grobner basis). For any /, g £ Rq define their S-pair as 



S(f,g) 



lcm(in> /, in> g) 



f- 



lcm(in> /, in> g) 



9- 



in> / 



Next define the normal form of an element /ei? with respect to a finite set of polyno- 
mials G = {gi, . . . , g r } £ Rq to be a polynomial 



for some unit u and polynomials such that in> / > in> NFg(/), in> / > in> a^j, and 
in>NF(j(/) is not divisible by any in>/j. Such a polynomial always exists and can be 
calculated explicitly using Mora's tangent cone algorithm [TT] (this is the local equivalent 
of the division algorithm). Buchberger's algorithm proceeds as follows: Starting with the 
generators of J, G — {fx, . . . , f s }, calculate NFc(5'(/j, fj)) for each pair i ^ j. If any are 
non-zero, add them to the set G and repeat the process, otherwise G is a standard basis 
for I. 

A standard basis G for I provides answers to many of the questions one might have 
about the local properties of J. For instance, / £ I if and only if NF G (/) = so a 
standard basis provides an algorithmic way to answer the ideal membership question. 
From a standard basis we can also calculate the Hilbert function of / which determines 
dimension of the component of J through b, and if b is an isolated solution it determines 
its multiplicity. 

Definition 3. The Hilbert function of ideal / with local order > is Hj : N — > N where 
Hj(d) counts the number of monomials with degree d that are not in in> /. 

Consider in> / in the lattice of monomials. For in > G = {mi, . . . , m r }, each monomial 
rrii cuts out the cone C mi of all its monomial multiples, and the monomials in in> I 
are exactly [J i C mi . The resulting picture is called a "staircase" (see Figure [[]) and the 
Hilbert function counts the monomials outside the staircase at each degree. Using the 
inclusion-exclusion principle and noting that C mi D C m . = Ci cm(mi)m .) we get an explicit 
combinatorial formula for the Hilbert function: 



fined this way, for fixed q the binomial coefficients are polynomial in p for all p > 0, 
and this polynomial has degree q. As a result, it is clear that the Hilbert function is 
described by a polynomial for sufficiently large degree. The regularity is bounded by 
deglcm(m!, . . . ,m r ) — n + 1, which is when all the binomial coefficients in the sum be- 
come polynomial. 

Definition 4. The g-corners of ideal / are the monomials that minimally generate in> /. 

The set of g-corners is uniquely determined by / and the local order >. The g-corners 
of / can easily be found from in> G for any standard basis G, and the g-corners completely 
determine Hj. 



NF G (/) 



uf ~ a\9i a r 9r 
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Figure 1: A zero-dimensional staircase, and a one-dimensional staircase. 

Claim 5. For sufficiently high degree, Hj can be expressed as a polynomial p. If p = 
then b is an isolated point of J. Otherwise the largest dimension component of J through 
b has dimension a + 1 where a is the degree of p. 

Unfortunately standard basis algorithms such as Buchberger are not numerically sta- 
ble. Small error in the calculation of b will cause large errors in the output. An alternate 
approach is to use the dual space instead, which can recover the same information as a 
standard basis and can be calculated with numerical stability. 



2 The Dual Space 

Considering R as a vector space, for each monomial x a there is a linear functional d° 
Ro — > k in i?o defined by 

Let Dq denote the vector space spanned by these monomial dual vectors. We will refer to 
D as the dual space of Ro even though it is technically a proper subset of Rq- By equipping 
Do with multiplication d a d 13 = d a+lS , it has a /c-algebra structure Do = k[di, . . . , d n ] where 
di is the dual element corresponding to x,. We give this ring a global monomial order >- 
which is the reverse of the order on R m , so if x a < x 13 then d a y d 13 . 

The dual space is sometimes defined in terms of differentials instead [12J. For / e Ro 

1 d ai d an 
ff*(f) = _ _ — f 



2=0 



Definition 6. The dual space of the ideal J, -Do [7]) is the subset of Do that annihilates 
/. The truncated dual space of J, D^[I], is the subset of D [I] of functionals with lead 
term of degree d or less. 

Theorem 7. Any monomial x a is in in> / if and only if the corresponding dual monomial 
d a is not in in^ D [I]. Equivalently, in^ D [I] = -Do[rn>/]- 
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Proof. For a proof in the case where / is zero-dimensional, see [9], Theorem 3.4. For 
positive dimensional I, fixing any a, let I' = I + m' a ' +1 where m is the maximal ideal 
(xi, . . . ,x n ). Note that V is zero-dimensional because Hr(c) = for c > \a\, so then 



x a G in> J' if and only if d a ^ in^ D [I']. Since x a has lower degree than any element of 
m M+i ; then x a G in> V if and only if x a G in> /. Additionally Dq [I] 



the elements of D 
and only if d' 



{d) 
o 

G in^ 



and m' a ' +1 have no terms in common. Therefore d° 
D [I). 



Dq[I'] because 
G m y D [I'] if 
□ 
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Figure 2: A standard basis (in blue) and dual basis (in black) for / = (x — y ,y ) 



Corollary 8. H^d) = dim Dj [/] - dim D^ 1 ^ [I] . 

Theorem 9. For any f G Rq, f G / if and only if p(f) = for all p G -Do[/|- 

Proof If / G I it is clear that p(f) = for all p G D [I]. Suppose f ^ I, and let 
G = {<7i, . . . , g r } be a standard basis of I. Then / can be expressed as 

/ = ttoig-i H h ua r g r + u NF G (/) 

where u is a unit, each a, is a polynomial, and in> NFc(/) ^ in> /. Let the lead monomial 
of NFg(/) be x a . Because u is a unit its lead monomial is 1, so the lead monomial of 
uNFcif) is also x a . By Theorem [7] there is p G D [I] with lead monomial d a . Due 
to the reverse nature of > and >-, p and uNFa(f) have only their lead monomial in 
common. Therefore p(aNF G (/)) ^ 0. Note that p{uaigi) = since ua^gi G / so p(/) = 
p(«NF G (/))^0. □ 

As a consequence, knowing a basis for the Dq [/] at each degree c up to some finite 
degree d provides some of the same information as a standard basis. In particular it 
reveals the values of the Hilbert function of / for all degrees up to degree d. Algorithms 
exist for finding a basis for the truncated dual space up to any particular degree. Two 
such algorithms are discussed below. Both reduce the problem to a system of linear 
constraints. Finding the kernel of a matrix can be done in a numerically stable way by 
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using singular value decomposition (SVD), which is what makes the dual space approach 
better suited for numerical situations. 

If I is known a priori to be zero- dimensional then Dq[I] has finite dimension, and it 
is possible to find an explicit basis for it. Finding the truncated dual at each successive 
degree, there will be some d for which D^[I] = D [I], and so dim D^ 1 ^ [I] = dimZ?Q [J], 
at which point we know the entire dual space has been found. The ideal membership 
question can then be answered for any element of Rq and the entire Hilbert function can 
be calculated (-ffj(c) = for all c > d). 

If / is not zero- dimensional, or if the dimension is not known, this strategy will not 
work. In general the dual space of I is not finite dimensional, so it is not possible 
to explicitly compute a basis for the entire dual space. The best we can do is find a 
truncated dual basis up to any finite degree d. The difficulty with this approach is that 
it's difficult to tell what degree d one needs to compute to in order to find all the relevant 
information about the ideal. Additionally, the truncated dual space can't be used for an 
ideal membership test. Given a polynomial /, even if p(/) = 0for ah>G D$[I\, it may 
still be that f £ I. 

Something that the truncated dual space can tell us is what the g-corners of / are up to 
degree d. Suppose all the g-corners of I are known up to degree d — 1. Then the g-corners 
at degree d are exactly the monomials missing from in^ [I] that are not multiples of 
the previously found g-corners. Computing [I] for successive d, if we could determine 
at what point all the g-corners of I had been found, then we could fully describe the 
Hilbert function Hj, since Hj is determined by the g-corners of /. We will present a way 
to do so. These methods will also provide a way to answer the ideal membership test for 
polynomials up to some fixed degree. 

2.1 Dayton- Zeng Algorithm 

A simple truncated dual space algorithm is one by Dayton and Zeng [I], using ideas of 
Macaulay [TU] . Given a finite generating set F, I considered as a vector space can be 
expressed as the span of the monomial multiples of the generators: 

/ = span{x a /| / G F, a G N n }. 

Proposition 10. Dq[I] is the set of junctionals p G satisfying p(x a f) = for all 
f G F and a G N n with \a\ + deg in> f < d. 

Higher degree multiples of the generators need not be considered when calculating 
D^[I] because they will not have any terms of degree d or less, so will always be or- 
thogonal to Dq [I]. Let Ad be the set of elements x a f as above. To find the subspace 
of Dq that annihilates all of Ad, we construct the Macaulay array M(F, d), which is 
the coefficient matrix of the elements of Ad- The matrix M(F, d) has entries in k with 
columns indexed by the monomials of Rq with degree < d and a row for each a G A d . If 
a is the ith element of A and x 13 is the j'th monomial with degree < d, then wiy = d l3 (a), 
with M(F, d) = (rriij). The truncated dual corresponds to the kernel of M(F, d). 
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Example 11. Let F = {x — y 3 ,x 2 } C k[x,y]. Then the Macaulay array M(F, 3) is 
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The kernel of this matrix has dimension 4, and a basis for it corresponds to the functionals 
1, d y , d 2 and d 3 + d x . These form a basis for the truncated dual space D^[I]. 



2.2 Mourrain Algorithm 

The second algorithm is due to Bernard Mourrain [12]. We define the "derivative" of a 
dual functional with respect to a given variable <9,. Let dj : D — > D be the linear map 
defined by 

[0 otherwise 

Note that djdj = djdj. Also for any / G Rq and p G Dq we have p{xif) = dip(f). It 
follows that if p G Dq[I] then d^p G Dq[I] for all i. In fact there is a stronger result: 

Theorem 12 ([12], Theorem 4.2). For any p G -Do? P £ 7) [7] */ an <^ on ^/ d?P G -Do[7] 
/or a// i and = /or a// f & F. 

The dual elements with lead term of degree d have derivatives which have lead term 
of degree d — 1 or less. This produces a way to build up D [I] degree by degree. Suppose 
0x,...,/3 r are a basis for Dq ^[1]. Then for p G D^[I], each derivative d^p can be 
expressed in terms of this basis so 

r 
3=1 

for some coefficients A* G k. It can be shown that 

r 

P = Y^ A )x x B j I n + X 2 x 2 0i I n H h A?x n & 

3=1 

so the elements of [71 are all linear combinations of the terms of the form xS~ I 
However not all linear combinations work. The fact that djd/p = d/d,p produces the re- 
lation 

r 

^(Ajd^- - Xjdtfj) = for each 1 < % < I < n. 

3=1 
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Since each dj/J,- is also in Dq [J], it can be uniquely expressed in the basis . . . , (3 r as 
dj/3j = Y7i=i A- Then the above equation can be broken down into the linear relations 

r 

for each 1 < i < I < n, and 1 < m < r. Finally p(fi) = for each generator /j produces 
another set of relations. We can build a matrix with a row for each of these constraints 
and columns corresponding to the coefficients A*. The kernel of this matrix corresponds 

to the space [I] . 

3 Homogeneous Ideals 

If / is homogeneous, then there is a criterion for deciding when all g-corners of / have 
been found when searching degree by degree. If / and g be homogeneous polynomials in 
Rq, then their S-pair is also homogeneous and 

deg S (/, g) = deg lcm(in> /, in> g) < deg / + deg g. 

In addition, if F is any set of homogeneous polynomials then NF p(S(f, g)) is also homoge- 
neous with the same degree as S(f, g). Suppose that F C / is a finite set of homogeneous 
polynomials with lead terms representing each of the g-corners up to degree d, but F is 
not a standard basis. By the Buchberger criterion, there is some pair f\, f'2 G F with 
g = NFp(S(fi, f'2)) 7^ 0. Then in> g is not divisible by any of the g-corners in in> F, and 
deg g < deg fi + deg fi < 2d, so there must be another g-corner with degree < 2d. 

Proposition 13. If I is a homogeneous ideal and C is the set of all g-corners of I up to 
degree d then either C is the set of all g-corners of I or there is an additional g-corner m 
with degm < max a 6eC - deg lcm (a, b) < 2d. 

So if finding bases for [I] at each c up to c = 2d reveals no g-corners with degree 
above d, then all g-corners have been found. We would like to extend this idea to the 
more general case of non-homogeneous ideals. Note that the bound of 2d can often be 
improved by taking max ai i> g c deg lcm(a, b) instead. 

Let R be the localization of k[t, Xi, . . . , x n ] by the maximal ideal (t, x 1 , . . . , x n ). For 
f E R let f h E R denote the homogenization of /. The ecart of / is the difference 
in total degree of the highest and lowest degree terms of /, or in other words the t- 
degree of the lead term of f h . For p e Dq let p h G k[dt,di, . . . ,d n ] = Dq denote the 
homogenization of p. Let ip : Dq — > Dq be the function that dehomogenizes with respect 
to t. So ip(q) = q\d t =i, and ip(p h ) = p. Abusing notation, we will also use ip to denote the 
dehomogenization function for i? . For I — (f 1: . . . , / s ), let / = (/f, . . . , f^) C Rq. Note 
that / is not the same as the homogenization of I and it will depend on the choice of 
generators of /. We fix a particular set of generators F = {f\, . . . , f s } from here forward. 
It is easy to see that ip(I) = I, regardless of the choice of generators. 

We take the local order on R to be some extension of the order > on R which is 
anti-graded and has t > Xi for every i. This ensures that for homogeneous g G Rq, 
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in> tp(g) = ■0(i n >fl')- The monomial order on Dq is taken to be the reverse of the order 
on R , so d t -< di for each i. Note that for homogeneous q G 7) this implies that 
myijj(q) = *jj(m y q). 

Theorem 14. If G = {g 1 ,...,g r } is a homogeneous standard basis of I, then ip{G) = 
{ip(g)i, ■ ■ ■ ,ip(g) r } is a standard basis of I . 

Proof. Any homogeneous g G 7 can be expressed as g — Yl,i a ifi-> so ^(flO = Z^VK *)/* ls 
in 7. Therefore ip(G) C 7. Moreover , 0( m > 9) = i n >' ^ / ; (fl , ) in> 7. For any / = X^i^/j 
7, let d be the maximum degree of all {bifi) h , and q be the integer such that t Ci (bifi) h 
has degree d. Then g = J2it Ci (bifi) h G 7 is a homogeneous degree d polynomial and has 
ip(g) = f and so ^(in > g) = in> /. Therefore ?/'(in > 7) = in> 7. For any m which is a 
g-corner of 7, t a m G in> 7 for some a. Taking a to be the minimum such value, t a m is 
a g-corner of 7. Therefore some g G G has t a m as its lead monomial, and m is the lead 
monomial of ip(g), so (in> ip(G)) = in> 7. □ 

Therefore we can find the g-corners of 7 by calculating 73^[7] for successive d, and 
using the stopping criterion for homogeneous ideals, and from this we can recover the 
g-corners of 7, which determines the Hilbert function Hj. 

Example 15. Let 7 be the ideal 

I = ( x 2 - xy 3 ,x 4 ) C C[x,y]( XtV ). 

All terms of the generators have degree 4 or less and the Hilbert function Hj(d) = 1 for 
d < 10. Finding the truncated dual of 7 at several degrees, one might be tempted to 
conclude that the Hilbert function stabilizes at 1. However, at d = 10 there is a new 
g-corner, and Hi(d) = 1 for all d > 10. A reduced standard basis of 7 is {x 2 — xy 3 ,xy 9 }. 
We look instead at the ideal 

I = (fx 2 - xy 3 , x A ) c C[f, x, y]( t ,x,y) 

which has reduced standard basis {t 2 x 2 — xy 3 ,x 4 ,x 3 y 3 ,x 2 y 6 ,xy 9 }. The g-corners of 7 
occur at closer intervals in degree. Beyond the highest degree of the generators of 7, no 
g-corner has degree more than twice that of the smaller degree g-corners. 

4 Eschewing Homogenization 

Although the method described in the previous section works, we would like to discover 
the g-corners of 7 without explicitly homogenizing the ideal, since this may introduce 
unnecessary numerical error to the process. Additionally, introducing an extra variable 
causes a significant increase in the computation time of the dual space algorithms, which 
we would like to avoid. To get around homogenization we can take advantage of the 
particular structure of 7. 

Definition 16. For / G 7? , the ecart of / is the difference in degree between the highest 
degree term and the lowest degree term of /. 
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Theorem 17. Let e be the maximum ecart of the generators fi, . . . , f s of I, and let q G Dq 
be homogeneous with lead term having d t -degree at least e. Then q G D [I] if and only if 
#7)6 A)[/]- 

Proof. Suppose ijj(q) G D [I]. For any / G /, q(t b f h ) = 4>{q){f) = if q and have 
the same total degree, otherwise q{t b f h ) = since they have no compatible monomials. 
Note that any element of / can be expressed as the sum of homogeneous polynomials of 
the form t b f h with / G /, and q annihilates all such polynomials, so q G -Do[^]- 

To prove the other direction, we use induction on the total degree d of q. For any d < e 
we have that ip(q) G D [I] implies q G D [I] vacuously since there are no functionals q 
with total degree d and lead term having <9 r degree at least e. Suppose for some d that for 
all homogeneous q G D [I] with degree at most d and lead term with d t - degree at least e 
that ip(q) G -D [/]. Fix some q G D Q [I] with degree d + 1 and lead term with ^-degree 
at least e. To show that ip(q) G D [I] it is sufficient to show i) that the first derivative 
diip(q) is in D [I] for each dual variable <9; and ii) that ip(<l)(fj) — for all generators fj. 

i) It is easy to check that differentiation with respect to any di commutes with deho- 
mogenization: diip(q) = ip(diq). The functional diq has total degree d and every term of q 
has c^-degree at least as large as the ^-degree of the lead term, which is > e, so the lead 
term of diq has ^-degree at least e, or d^q = 0. Therefore diq G D [I] so di(ip(q)) G D [I]. 

ii) ip(q)(fj) = d^q(t b fj') for any values of a and b for which the total degrees of 
<9"g and t b f^ are equal. If a > b then &lq{t b f^) = since every term of d?q has dt- 
degree at least a + e and every term of t b f^ has i-degree at most b + e. If a < 6 then 

= q(t b ~ a f?) = since g G D [/| and t 6 " // e/- □ 

Corollary 18. Lei e 6e t/ie maximum ecart of the generators of I. Every g-corner of I 
has t-degree < e. 

Proof. If t a x a is a g-corner of / then x a G in> I by Theorem [TH and so d a ^ in^ D [I}. 
Therefore d b d a ^ in^Dol-T] for all b > e, so G in> / for all b > e. Since a is the 
minimum value for which t a x a G in> J, it must be that a < e. □ 

At and above t-degree e, the dual space of I looks just like the dual space of I (after 
dehomogenizing) . At and below t-degree e is where g-corners of / may occur and this 
information will be used to decide what degree to calculate the dual space up to. Let 
Dq[I] denote the subspace of D [I] with degree exactly d. 

Corollary 19. Dq[I] C ^(Dq[I\) and the subspace of i/j(Dq[I]) of elements with lead 
term of degree d — e or less is equal to Dq e ^ [I] . 

Proof. For any p G [I] let q be the homogenization of p to degree d, that is q = d^p h 
with a chosen so that q has degree d. For any g G J, q(g) = trivially if g does not have 
degree d, and otherwise q(g) = p(ip(g)) = since ip(g) G /. Therefore q G -Do[/], which 
proves the first part of the statement. An element in i/j(Dq[I]) with lead term of degree 
< d — e is the dehomogenization of an element p G Dq [I] with lead term with <9f-degree 
> e, so tp(p) G D [I). □ 
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Figure 3: Here I = (x — y 3 , x 3 ). On the left is the monomial lattice for J, sliced at d = 6. 
All g-corners occur at or below the plane of t-degree 2. On the right is the slice after 
dehomogenization where black dots represent lead monomials of i/j(Dq[I]) 



5 The Sylvester Dual 

Definition 20. The Sylvester array S(F, d) for a set of generators F of ideal I, and a 
degree d is the coefficient matrix of all monomial multiples of the generators in F, x a f 
such that every term of x a f has degree d or less. The columns correspond to each of the 
monomials up to degree d. 

The Sylvester array S(F, d) is similar to the Macaulay array M(F, d) but instead of 
having a row for every monomial multiple of a generator that has any terms of degree 
< d, it only includes the ones that have all terms of degree < d. The kernel of M(F, d) 



corresponded to D^[I]. The kernel of S(F, d) also defines a subspace of Dq 1 ' , which will 



(d) 



be denoted Sq [F]. Note that unlike the truncated dual space, Sq^F] depends on the 
set of generators for I. Also unlike the truncated dual space, it is not generally true that 



(<*)[ 



S^[F]CS^'[F}. 

Example 21. As in the Macaulay array example (Example ITT]) let F = {x — y 2 ,x 2 } C 
k[x,y}. Then the Sylvester array S(F, 3) is 



(d+l) r 



1 Orr 8„ 



x-y 

x 2 

x(x 2 
y(x 2 



3 /0 





1 














d 2 

X 



1 






d y d x 







d 2 

v 







dl d„dl d 2 d x 





1 








1 










a 3 

v 

-1\ 





The kernel of this matrix, Sq [F], has dimension 6, with basis {1, d y , d y d x , d 2 , d 2 d x , d 3 + 
d x }. Note that the rows of S(F,d) are a subset of the rows of M(F, d), so DQ d \l] is 
contained in Sq[F\. 
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Theorem 22. S ( d) [F] = i>{D$[I]). 

Proof. Dq[I] is exactly the set of degree-d homogeneous functionals in Do that annihilate 
the degree-d homogeneous polynomials of /. The space of all degree-ci homogeneous 
elements of / is spanned by the elements of the form t a x a f h where / G F is part of the 
generating set of I and a + \a\ + deg f h = d. Note that deg f h is equal to the maximum 
degree of any term in /. Dehomogenizing everything, iP(Dq[I]) is the set of functionals in 
Dq that annihilate all polynomials of the form x a f such that all terms have degree < d. 
These polynomials are exactly the ones which form the rows of S(F, d). □ 

This relationship provides an alternate way to prove Theorem [171 

(Alternate proof of Theorem\Tl\) . Let q G D be homogeneous with total degree d and 
lead term with <9 t -degree > e. This implies that if>(q) has all terms of degree d — e or 
less. Suppose q G Dq[I], so then ifj(q) G Sq[F] by the previous theorem. This means 
if)(q) annihilates each x a f with all terms of degree d or less and / G F. For x^f with 
some term of degree > d, the degree of the lead term of ar / must be > d — e because 
the ecart of / is at most e. Therefore if>(q) also annihilates x@f, since they have no terms 
in common. ip(q) annihilates all terms of the form x a f for any x a and any / G F, so 
1>(q)ED [l\. 

Suppose if>(q) G D [I]. Then ip(q) annihilates all x a f for / G F, which implies 
^(q)eS^[I\, so qeD [I\. □ 

[F] can be calculated at each degree d without homogenizing, and captures all the 
information of the homogenized dual space. In the lattice of monomials, in^ Z)q [J] can 
be considered as a slice of in>_ Dq[I] at degree d. For each g-corner m of I, the monomial 
multiples of m will be missing from in>_ -Dq[/]. This slice of the cone generated by m will 
appear as a truncated cone of missing monomials in Sq^[F], starting at the monomial 
if)(m) and extending out to all multiples of if>(m) up to degree d — a where a is the t- 
degree of m. The monomials missing fromS^F] are the union of all the truncated cones 
generated by the all the g-corners of I up to degree d. 

In Dq [I] there is also a cone of missing monomials at if) (m) for each g-corner m of J, 
but in this case the cone extends all the way to degree d. 

Suppose {gi, . . . ,g r } is the standard basis of /. The subspace of I with degree d is 
spanned by polynomials of the form mgi where m is any monomial with degree d — deg g^. 
For any particular g iy the possible values of if){m) are all the monomials in R up to degree 
d — deg g^. In the lattice of monomials of Rq, the possible values of in > ('0(m<7j)) form a 
truncated cone, starting at in> (^(^i)) an d extending out to all multiples up to degree 
d — (deg in > ((7 i ) — deg in>(^(pj))). The value of degin>((7j) — deg in>(^(^j)) is the t-degree 
of in>((7j) and it is at most e. The monomials excluded from in>(?/>(Z)g[/])) are the union 
of the truncated cone from each 'm > (if}(gi)). In contrast the monomials excluded from 
iiiylD^lI]) are all multiples of in>(^i(gri)), so this picture is similar but the excluded 
cones extend all the way to degree d. 

If all g-corners of I are known up to degree d — 1, then it can be calculated exactly 
which monomials should be missing in in^ S { d) [F] if there are no additional g-corners at 
degree d. By calculating a basis for S^fF], whichever additional monomials are missing 
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Figure 4: Black dots represent lead monomials in Sq [F] for d — 5 and d = 6 with 
.F = {x 2 — xy 3 ,x 4 }. The missing monomials form truncated cones (blue triangles). At 
d = 6 a new g-corner is discovered at x 3 y 3 . 



must be new g-corners of I with degree d. Therefore by calculating in^ Sq [F] for each 
successive d it is possible to discover the g-corners of I at each degree. We can determine 
when all the g-corners of / have been found using Proposition [13] The g-corners of / 
determine the g-corners of /, which fully determines the Hilbert function Hi. The process 
of finding a basis for 5"q^[F] also produces the truncated dual space of /, since at each d 
the set the elements of [F] with lead term at most d — e is ^ [I] . 

Algorithm 23. 

Inputs: generators F = {fi, . . . , f s } of ideal /. 

Outputs: monomials C = {ci, . . . , c r } corresponding to g-corners of /. 
B := {}; List of pairs of a g-corner and a degree. 
d := 0; 

dmax ■= 2 maxi{deg in> fi}; 
while d < d max do 

build Sylvester array S(F, d); 

Sdual := basis for ker S(F, d); 

reduce Sdual so each element has unique lead term; 
B d := {}; 

for monomials m ^ in > Sdual with deg m < d do 

if for all (q, dj) G B, either q { m or deg(m/c,) > d — di then 

append (m, d) to 
end if 

end for 

if B d ^ {} and d max < 2d then 
end if 

append B d to B; 
d 4- d + 1; 
end while 
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C := list of monomials from pairs in B; 
return C; 

Given the output C of this algorithm, to find the g-corners of /, simply remove all 
monomials q G C that are divisible by some other monomial in C. 

Remark 24. When dealing with numerically obtained inputs, the algorithm should use 
singular value decomposition to calculate a basis for the kernel of the Sylvester array. 
This method is numerically stable, while using Gaussian elimination is not. 

Sq [F] reveals not only the g-corners of J at degree d (and therefore the g-corners 
of I), but it can also be used to find the corresponding standard basis elements of /. 
Additionally S^[F] can answer the ideal membership test for polynomials with all terms 
of degree d or less. These two facts are encompassed by the following proposition. 

Proposition 25. If f G Ro is a polynomial with no terms exceeding degree d andp(f) = 
for all p G Sq[F], then f G I. For each monomial d a of degree < d that is not in 
in^ Sq [F] , there is some polynomial f G Ro satisfying the above conditions with in> / = 
x a . 

Proof. Suppose / G -Ro is a polynomial with no terms exceeding degree d and p(f) = 
for all p G S$ [F]. Let g be the homogenization of / to degree d (i.e. g = t a f h where 
a = d — deg f h ). Then q(g) = for all q G Dq[I] so g G I. Therefore / = ip(g) G J. 

If d a with degree < d is not in in^ [F] , then its homogenization d^d a to degree d 
is not in in^ D Q [I] so there is some homogeneous g G / with in> g = t b x a . Therefore ip{g) 
has lead term x a and is annihilated by Sq [F] . □ 

Supposing a basis for Sq^ [F] has been calculated, build the coefficient matrix of these 
basis elements with columns for each of the monomials up to degree d. The kernel of 
this matrix corresponds to the polynomials in R with all terms of degree < d that 
are annihilated by Sq [F]. Let mi, . . . ,m r be the set of g-corners of /. If monomial 
r/ij G -Do has degree d, there must be some polynomial hi in this kernel with in> gi = 
4>{hi). Collecting the polynomials found this way for each g-corner of I produces a set 
H = {hi, . . . , h r } C /, with in> H = in> / so H is a standard basis of /. 

Algorithm 26. 

Inputs: Basis for the Sylvester dual Sdual at some degree d, and a g-corner c of I found 
at degree d. 

Outputs: Polynomial p G / with in> p = c. 

monomials := list of monomials m G Ro with degm < d and m < c; 

M := coefficient matrix of elements in Sdual, with columns only for the monomials in 

monomials; 

K := kerM; 

return an element p of K with in > p = c; Such an element is guaranteed to exist. 

Example 27. We continue Example [15] with F = {x 2 — xy 3 ,x 4 } and / = (F), and run 
through the algorithm for finding the g-corners of I and a standard basis, this time using 
Sylvester arrays instead of homogenizing the generators. For d = 0,1,2,3 the Sylvester 
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arrays are empty because there are no multiples of x 2 — xy 3 or x 4 which have all terms of 
degree 3 or less. Therefore bases for the first four Sylvester dual spaces are 

S^iF]: {!}, 

{l,d y ,d x }, 
S^[F}: {l,d y ,d x ,d 2 ,d y d x ,d 2 x }, 
S ( V[F] : {l,d y ,d x ,d 2 y ,d y d x ,d 2 x ,d 3 ,d 2 y d x ,d y d 2 x ,d 3 x }. 

At degree 4, the Sylvester array S(F, 4) has rows for x 2 — xy 3 and x 4 . The kernel of 
this matrix [F] has basis 

S^IF] : {1, d x , d y , d y d x , d 2 , d 3 x , d y d 2 x , d 2 d x , d 3 , 

d y d 3 x) d 2 y d 2 x) dld x + d 2 x) dl}. 

The monomials d 2 and d 4 are both missing from the set of lead monomials of these 
basis elements. Since there were no previous g-corners found, each of x 2 and x 4 must be 
the dehomogenization of a g-corner of /. They are recorded as potential g-corners of / 
along with the degree they were found at, which is 4. To find standard basis elements 
corresponding to these g-corners, we construct the coefficient matrix of the basis for 
.Sq 4 '' [F], and try to find elements of the kernel with lead terms x 2 and x 4 . The polynomials 
x 4 and x 2 — xy 3 may be produced this way. 

At degree 5, S(F, 5) has rows for x 2 — xy 3 , x 4 , x 3 — x 2 y 3 , x b , x 2 y — xy 4 , x 4 y. A basis 
for the Sylvester dual is 

S^[F] : {l,d x ,d y ,d y d x ,d 2 y ,d 2 y d x ,d 3 ,d y d 3 ,d 2 d 2 x , 

d 3 d x + d 2 x , d 4 , d 2 d 3 , d^d 2 , d 4 d x + d y d 2 , dl). 

Since x 2 and x 4 were g-corners found at degree 4, the multiples up to 1 degree higher will 
be missing at degree 5. This accounts for all the monomials missing from this basis for 
S^[F], which are d 2 , d 3 , d y d 2 , d 4 , d x and d y d 4 , so there are no new g-corners here. 
At degree 6, a basis for the Sylvester dual is 

S^[F] : {1, d x , d y , d y d x , d 2 , d 2 d x , d 3 , d 3 d x + d 2 , 
d 4 y ,d 2 y d 3 x ,d 3 y d 2 x + d 3 x ,d 4 d x + d y d 2 x , 

d 5 y ,d 4 y d 2 x + d y d 3 x ,dp x + d 2 d 2 x ,d e y }. 

The monomials missing from the lead terms are d 2 , 8%, d y d 2 , d^, d y d 3 , d 2 d 2 , d x , d y d 4 , 
d x , d y d x , d 2 d 4 and d^dl. All but d^dl corresponds to a g-corner x 2 or x 4 or a multiple 
of these by a monomial up to degree 2. Therefore x 3 y 3 is a new g-corner, which we store 
along with the degree 6 at which it was found. (See Figure 4 for a diagram of this step.) 
A corresponding standard basis element is x 3 y 3 . 

Continuing this process at each degree, the g-corner x 2 y & is found at d — 8 and 
xy 9 is found at d = 10. Corresponding standard basis elements are x 2 y 6 and xy 9 . No 
additional g-corners are found searching up to degree 20, so all the g-corners of / are 
among x 2 , x 4 , x 3 y 3 , x 2 y 6 , xy 9 . The monomials that are multiples of other monomials in 
the list can be dropped, leaving x 2 and xy 9 . The corresponding standard basis elements 
found for these g-corners were x 2 — xy 3 and xy 9 so this is the reduced standard basis for 
J. Finally the Hilbert function can easily be recovered from the set of g-corners, which is 
ifj(d) = 2 for 1 < d < 10, and H T (d) = 1 for d = and d > 10. 
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We have implemented Algorithm [23] for finding the g-corners of an ideal using the 
Sylvester dual, and this algorithm for recovering a standard basis of the ideal, in the 
computer algebra system Macaulay2. This implementation is contained in the package 
"NumericalHilbert," which can be found at 
http : / / people . math . gatech . edu/~ rkrone3/NHcode . html 

[F] can also be calculated using a variation of the algorithm presented by Mourrain. 
The Mourrain algorithm works by finding the elements of D whose derivatives are in the 
span of the previously found dual basis elements, and which annihilate the generators 
fi, . . . , fg. Here we can take advantage of the fact that there exists a homogeneous basis 
of the truncated dual space. Each of the derivatives of a degree d dual element has degree 
d — 1, so to find the degree d elements we need only to consider the elements found in the 
preceding step. The Macaulay2 package also contains an implementation of this version 
of the algorithm. Using the Mourrain algorithm is more efficient than using the Sylvester 
array strategy in most cases because the matrices involved grow with the dimension of 
the dual space rather than with the dimension of the entire ring up to the given degree. 

The algorithm produces: 

• a basis for the dual space truncated to the degree that the algorithm needed to 
calculate up to 

• the g-corners of the ideal (or optionally a standard basis of the ideal) 

• a bound on the regularity of the Hilbert function 

• the values of the Hilbert function up to the regularity bound 

• the Hilbert polynomial which defines the values above the regularity 

Note that the exact regularity of the Hilbert function can easily be computed from this 
data by comparing the returned values of the Hilbert function to the Hilbert polynomial 
at each degree below the bound. 

Example 28. In this example we demonstrate the functionality of the Macaulay2 pack- 
age. Let / be the ideal defined by the system Cyclic4, given by the following polynomials 

in C[:ci,X2,X3,:c 4 ]: 

F = {x\ + x 2 + x 3 + x 4 , Xix 2 + x 2 x 3 + x 3 x 4 + £4X1, 

X1X2X3 + X1X2X4 + £1X3X4 + X2X3X4, X1X2X3X4 — 1}. 

The variety consists of two irreducible curves, along with 8 embedded O-dimensional com- 
ponents. One of the embedded components is (—1, 1, 1, —1). However suppose we were 
not aware of this, and only had an approximate numerical value for this point. Below we 
use an approximation with small error that was obtained using a numerical solver. 

11 : loadPackage "NumericalHilbert"; 

12 : R = CC[x_l . .x_4, MonomialOrder=>{Weights=>{-l,-l,-l,-l}},Global=>false] ; 

13 : F = matrix{{x_l + x_2 + x_3 + x_4, 

x_l*x_2 + x_2*x_3 + x_3*x_4 + x_4*x_l, 

x_2*x_3*x_4 + x_l*x_3*x_4 + x_l*x_2*x_4 + x_l*x_2*x_3, 
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x_l*x_2*x_3*x_4 - 1}}; 

14 : P = {-1.0-.53734e-17*ii, 1 . 0- . 20045e-16*ii , 

1.0+.89149e-17*ii, -1 . 0+ . 18026e-17*ii} ; 

15 : dualInfo(F, Point=>P, Tolerance=>le-4) 

o5 = ({1, -lx +x,-lx + x , . . . }, 
13 2 4 

2 

{x , x , x , x x }, 5, {1, 2, 1, 1, 1}, 1) 

1 2 3 3 4 

The output at o5 is a sequence consisting of the information listed above the start of 
the example. First is a basis for [I] but we do not reproduce the full output here in 
the paper for reasons of space. Shown is only D^[I\. Note that the dual elements are 
written in terms of the original variables of Ro even though this is an abuse of notation. 
Next in the sequence is the list of g-corners generating in> I, x\,X2-,x\,X2X±. The last 
two entries of the sequence are the first five values of the Hilbert function (1,2,1,1,1) 
followed by the Hilbert polynomial, which is 1. We see from the output that the Hilbert 
function is Hj(d) = 1 for all d except -ffj(l) = 2. The fact that the Hilbert polynomial 
is a non-zero constant indicates that indeed the point in question sits on a 1-dimensional 
component of the variety. All of these results agree with the values that are obtained by 
symbolic computation at the point (—1,1, 1,-1), but we did not need to know the exact 
value of this point. 

Example 29. We give another example, this time where the point of interest is not 
rational so it must be approximated. Define I to be the ideal generated by F = {{xf + 
x l + x i ~~ iX^i — x 2 ), (xi — X2) 3 } C C[xi, X2, X3]. The variety is the plane x\ — X2 = 0, 
but there is an embedded curve which passes through the point (\/2, \/2, 0). We will use 
a numerical approximation of this point to find the Hilbert function of the ideal localized 
here. 

11 : loadPackage "NumericalHilbert" ; 

12 : R = CC[x_l..x_3, MonomialOrder=>{Weights=>{-l , -1 , -1}} , Global=>f alse] ; 

13 : F = matrix{{(x_l~2 + x_2"2 + x_3~2 - l)*(x_l - x_2) , 

(x_l - x_2)~3}}; 

14 : P = {0.7071068, 0.7071068, 0}; 

15 : dualInfo(F, Point=>P, Tolerance=>le-4) 

o5 = ({1 , x , x , lx , . . . }, 
12 3 

2 2 

{x , x x }, 4, {1, 3, 5, 6}, i + 3) 
1 12 

We again truncate here the list of basis elements of the dual space of the ideal for 
space reasons. The algorithm finds that 4 is a bound on the regularity of the Hilbert 
function, and the values of Hj(i) for % — 0, 1,2,3 are 1,3,5,6 respectively. Beyond that 
point, Hj(i) = i + 3. The fact that the Hilbert polynomial is linear indicates that the 
largest component passing through (\/2, v^2, 0) has dimension 2. 
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